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In cells, helicase translocation along nucleic acid is essential for many biological processes. How- 
ever, so far, the mechanism of this translocation is not fully understood. Recent studies show that 
helicase might translocate through two processes, active process and passive process, with different 
translocation rate. In this study, a mechanochemical model including such two processes is pre- 
sented for ribosome translocation along messenger RNA during translation. In which, each of the 
two processes consists of two sub-processes, translocation factor binding sub-process and the follow- 
ing mechanochemical sub-process in which ribosome makes a forward translocation step. Ribosome 



' switches stochastically between these two processes with external force dependent rates. By this 

model, we found that, with the increase of external force, the mean translocation rate of ribosome 
fT^ ' increases from one lower limit to another upper limit, and both of these two limits increase with 

■ concentrations of the translocation factors. Under high external force, ribosome translocates mainly 

through active process. At saturating concentration of translocation factors, the translocation is 
limited by mechanochemical sub-processes. 
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I. INTRODUCTION 



^ . In many biological events, such as replication, recombination and repair, double-stranded nucleic acid, i.e. DNA 
\^ • or RNA, should be firstly unwound to provide single-stranded template [lj-|6|- In this process, using the energy 
' released from ATP hydrolysis helicase translocates along one strand of the double helix to unwind the nucleic acid 



[6 - 12 1 . Single molecule experiments found that the translocation rate of helicase along nucleic acid depends on the 



stability of base pairs (GC or AU base pair), the external force (which can be applied by optical tweezers assay 



vn 
o 

or flow stretching assay), and the concentrations of translocation factors, such as EF-G, EF-Tu and ATP |13l42l| 

The translocation of helicase along nucleic acid might be active (in which the helicase acts as a strong molecular 

motor which converts chemical energy from ATP hydrolysis to unzip a duplexed template with high efficiency) or 

^ passive (in which helicase o nly can translate to the next base after the thermal fluctuation induced opening of the 
Cq — — . . . . . . 



adjacent base 
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which helicase o nly 
,se pair) 
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. So far, many models have been established to describe this unwinding process 



In this study, similar as in the description of dynamics of microtubules |35l.l36l|. a two-process model will be presented 



for ribosome translocation along messenger RNA (mRNA) during translation. In which, ribosome translocates along 
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mRNA by either active process or passive process. During translocation, ribosome switches stochastically between 
these two processes, and the switch rate from passive process to active process depends on external force. Each process 
consists of two sub-processes: the translocation factors (EF-G»GTP and EF-Tu»GTP»aa-tRNA) binding sub-process, 
and mechanochemical sub-process during which ribosome will make a forward translocation step (i.e. one base pair). 
This model fits well to the recent experimental data obtained in [23|. Based on this model, properties of the ribosome 
translocation along mRNA, including its mean translocation rate and mean dwell time in one translocation cycle, are 
detailed discussed. 



II. MODEL 



As discussed in recent references 



m 



23|, the translocation of helicase along nucleic acid (DNA or RNA) 



might through two processes: active process and passive process, see Fig. [TJ The main difference between these 
two processes is that, in active process, after completion of the necessary translocation factor binding processes, the 
junction of the two strains of the double-strained nucleic acid is effectively destabilized, and so the energy barrier 
for helicase to translocate to the next base is greatly reduced. But in the passive process, the helicase has to 
overcome a significant energy barrier, and so more time will be spent in waiting for the junction opening by thermal 
fiuctuation. Biochemically, the free energy AG needed to open one base pair consists of the base-paring energy AGbp, 
the destabilization energy of enzyme AGd and template stretching force AGp, AG = AGbp — AGd — AGf. 

In both active process and passive process of ribosome translocation along mRNA, translocation factors EF-G»GTP 
and EF-Tu»GTP»aa-tRNA are needed 23, 3^. Therefore, there are at least two sub- processes in each of the two 
processes. For simplicity, we assume there are only two sub-processes in each translocation cycle: one translocation 
factor binding sub-process, and one mechanochemical sub-process during which one spatial mechanical translocation 
step is completed. Generally, all the four sub-processes (see Fig. [T]): 1 — > 2, 2 ^ 1 and 1' ^ 3, 3 — > 1' might 
be reversible, for simplicity we assume in our model that they are all irreversible, with rate constants ki , fc2 and 
k'ljk^ respectively. 1 — > 2 ^ 1 represents one active translocation step, and 1' — s> 3 — > 1' represents one passive 
translocation step. As discussed in ^3], 2 — > 1,3 — !> 1' can be regarded as ribosome translocation through an open 
and closed nucleic acid fork, respectively. So kg, depends only on the junction base pair, and ^2, kg are independent of 
external force. To reduce model parameters (and based on numerical tests), we assume that only the switch rate from 
passive process to active process is external force dependent, and the rates of translocation factor binding sub-process 
in both active (1 2) and passive (1' — ?> 3) processes are equal to each other, ki — k[, which depend only on the 
concentrations of translocation factors. Generally, the binding of translocation factors to ribosome might generate 
conformational change of the corresponding domain, which then might cause a little mechanical translocation. Here, 
for simplicity and to reduce the model parameters, we assume that the corresponding rates ki,k[ are independent 
of external force. Which means the whole mechanical step is made in the following mechanochemical sub-process. 
Meanwhile, we assume the switch rates between states 1 and 1' are the same as those between states 2 and 3, i.e. 
fcg — ko and k'^ — kc- Under all these assumptions, the model illustrated in Fig. [T]is then reduced to the simple one 
illustrated in Fig. [D^a). To some extent, the model depicted in Fig. [UJa) is equivalent to the one depicted in Fig. 
O^b). In this study, we will only discuss the model depicted in Fig. HJa), since more information can be obtained from 
it but with the same parameter number as the one depicted in[2Jb). 
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III. RESULTS 



Mean translocation rate. 



Let pi,P2 and pi,p2 be the probabilities that ribosome in states 1, 2 and 1', 3 respectively, then one can easily 
show that they are governed the following master equations 

dpi/dt =fcoPi + k2P2 - (fci + kc)pi, 
dp2/dt ^koP2 + kipi - {k2 + kc)p2, 
dpi/dt =kcPi + k3p2 - {ki + ko)pi, 



dp2/dt =kcP2 + kipi - (fca + ko)p 



2- 



At steady state, one can easily get 

Pi — i^o 

P2 = koki{ki + ks + kc + ko)/A, 

(2) 

Pi - kc[k2{k3 + ko) + ksiki + /cc)]/A, 
P2 = kcki{ki + k2 + kc + ko)/A, 

with the normalization constant A = (fco + kc)[{ki + k2){ki + ^3) + kc{ki + k^) + ko{ki + fc2)]- The mean translocation 
rate of ribosome can be obtained as follows 

J = ki{pi + pi) ^ k2P2 + k3P2- (3) 

The ratio of probability p := pi + P2 that ribosome stays in active process to probability p '■= pi + P2 that in passive 



process is p/p = ko/kc- There is evidence that the ratio of switch rates ko/kc increases with external force 23[, so 
under high external force, ribosome will translocate mainly by active process. One can easily show that the mean 
translocation rate of ribosome in active process is Ja = + ^2)7 and the mean translocation rate in passive 

process is Jt — kikj,/{ki + fc3). For low concentrations of translocation factors, i.e. ki is small, the mean translocation 
rate J « fci, so the translocation is limited by the translocation factor binding process 1 — > 2 or 1' — s> 3. But for 
high concentration cases, i.e. ki is large enough, J k, (fc2fco + k'ikc)/{ko + fcc), so the translocation is limited by 
the following mechanochemical process 2 — >■ 1 or 3 — >■ 1'. One can also show that, the ratio {pi + pi)/ij>2 + P2) 
decreases with fci, which means that for low concentration of translocation factor cases, ribosome mainly stays in 
states 1 or 1', but for high concentration cases, it is mainly in sates 2 or 3. In one word, the external force determines 
which process, the active one 1 ^ 2 ^ 1 or the passive one 1' ^ 3 ^ 1', is the dominant process, while the 
concentrations of translocation factors determine which sub-process, the binding sub-process 1 — > 2 (1' — > 3) or the 
following mechanochemical sub-process 2 — > 1 (3^1'), limits the whole translocation of ribosome. 



B. Mean dwell time. 



The mean time that ribosome spent in one translocation cycle can be obtained as follows 

T = piTi + P2T2 + PlTy + P2T3. 



(4) 
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Where Ti,T2,Ti/,T3 are the mean times spent by ribosome to complete one translocation cycle with initial states 1, 
2, 1', 3, respectively. Ti satisfies the following equations {36, 

1 fcl 



38 



Tiv = 

1^13 = 



fcl + fcc fcl + fc, 
1 fco 



fcl + fc, 

fcl 

fcl + fco ' fcl + fco ^ ' fcl + fco 

1 fc, 



fc2 + fcc fc2 + fc, 
1 , fco 



■T12 

■ri + 

■Tl3, 
-^12. 



■Til', 



13, 



(5) 



fc3 + fco fc3 + fcc 

Where Tu is the mean first passage time of ribosome from state i to state 1 or 1' of the next translocation cycle [see 
Fig. (Ha)]. One can show that Ti = Si/S, with Ei = fci(fci + fc2 + fcc + ko){k^ + fcc + fco) + (fcc + fco)(fc2fc3 + fc3fcc + fc2fco) 
and E — fci(fc2fc3 + fcsfcc + fc2fco)(fci + fcc + fco)- The expressions of T2,ri',T3 can be obtained similarly. Ty — E2/E 
with E2 = fci(fci + fcs + fcc + fco)(fc2 + fcc + fco) + (fcc + ko){k2k^ + fc2fco + fccfcs)- 

(fcl + fc3)fcc + (fcl + fc2)(fc3 + fco) 



and 



To 



T3 



fcl(fc2fc3 + fcsfcc + fc2fco) 

(fcl + fc2)fco + (fcl + fc3)(fc2 + fcc) 
fcl(fc2fc3 + fc2fco + fc3fcc) 



(6) 



(7) 



C. Ribosome translocation along mRNA. 



In the following, we will use the above model to discuss the external force and translocation factor concentration 
dependent properties of ribosome translocation along mRNA, which has been recently studied experimentally in [2^. 
As mentioned before, the switch rates of ribosome between active process and passive process depend on energy barrier 
of the strains junction AG = AGbp — AGd — AGp- Under the assumption that AGbp and AGd are invariable in 
the translocation, the switch rates satisfy fco/fcc — fc^/fcc exp((5AGi?/fcBT). Here < (5 < 1 is an energy distribution 
factor, and numerical calculations indicate i5 < 1. For simplicity and reduction of model parameters (and also due to 
the test of numerical calculations), we assume fcc = fcj? is independent of external force. So 

fco = fc°exp(,5AGF/fcBT). (8) 

Another method which is usuall y u sed to account for the external force F dependence of transition rate is fc = 

This will be used in this study, i.e, we assume the external force 



36, 



4C 



fcoexp(i^/Fo) with Fq a constant 
dependence of switch rate fco is as follows. 
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ko^k°eMF/Fo). (9) 

In fact, the methods ([5]) and © are equivalent to each other for ribosome translocation process studied here. The 
reason is as follows. From the force-extension relation of DNA measured in 42] [or see|3l^a)], the energy AGf can be 
well approximated by 16F — 8 for F < 60 pN, see Fig. ^h) in which the thin solid line is from AGp = 16i^ — 8 and 
the thick dotted line is from AGp = z{F)dF with z{F) the extension as plotted in Fig. EJa). So, by expressions 
(|8]) and ([9]), one can easily show that 



X ^bT ^0 ,0 1 

= , fc" — fc" exp 

I6F0' ° ° \2Fq 



(10) 
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In the discussion below, we will always use expression and the same results can also be obtained by method ([5]). 

Generally, the rate of translocation in passive process should depend on the type of base pairs of the strains 
junction. Again, for the sake of simplicity, we always keep it constant [2^. Similarly, the transition rate ^2 in 
active process is also assume to be a constant. As studied in [23^, there are at least two binding processes during 
sub-processes 1 -t- 2 or 1' — > 3, EF-G»GTP binding and EF-Tu»GTP»aa-tRNA binding. For simplicity, we assume 
there are only these two translocation factor binding processes (or all other translocation factors are held at saturating 
concentrations), and both of them are irreversible. Let fcn = fc]']^[EF-G] and ki2 = fcj'jiEF-T] be the binding rates of 
EF-G»GTP and EF-Tu»GTP»aa-tRNA respectively. One can easily show that 

fcnfci2 _ fc°ifcP2[EF-G][EF-T] 
kn+ki2 A:0i[EF-G] + fc?2[EF-T]- ^ ' 

With parameter values listed in Tab. HI the theoretical results of ribosome translocation rate J versus external 
force F and concentrations [EF-G], [EF-T] are plotted in Fig. 21 where the data points are from [23] for hpValGcso 
mRNA. These plots indicate our model can explain the experimental data reasonably well. In Fig. [SJa), the curves 
for translocation rate J versus force F at saturating [EF-T] but various values of [EF-G] are plotted, and in Fig. \^c), 
the curves for translocation rate J versus [EF-G] at saturating [EF-T] but various values of force F are plotted. From 
these two figures, one can easily see that, the translocation rate J increases monotonically with both force F and 
concentration [EF-G], but with one upper limit and one lower limit. Similar results can also be found for different 
values of [EF-T], see Fig. [5jb)(d). As mentioned before, at saturating concentration of translocation factors, the 
translocation of ribosome is limited by the mechanochemical sub-process 2 — > 1 (or 3 — > 1'), but at low concentration, 
it is limited by the translocation factor binding sub-process 1^2 (or 1' 3). Meanwhile, under high external 
force, ribosome translocates mainly through active process, but under low force, it is mainly through passive process. 
Actually, this can be demonstrated by Fig. [HI From Fig. [D^c), one sees that for any external forces, the probability 
Pi+Pi of ribosome staying at state 1 or 1' decreases with concentration [EF-G], which is because that the translocation 
factor binding rate fci increases with [EF-G]. However, the curves in Fig. |6l^a) indicate that pi + pi increases with 
external force F, this is because that, with the increase of F ribosome will more like to reach state 1 (or 1') by 
sub- process 2 — > 1, which is fast than sub-process 3 — !• 1' (see Tab. [Jfor values of k2 and h^). This also can be seen 
from Fig. [Hb) that, the probability pi +p2 that ribosome translocates through active process increases monotonically 
with force F. Actually, numerical calculations show that this increase is almost exponential. The plots in Fig. IH^d) 
indicate pi -^-p2 is independent of concentration [EF-G] . In all the plots of Fig. |6l the concentration of [EF-T] is always 
held at a saturating concentration. Similar results can be obtained if [EF-T] is varied. Since all the sub-processes in 
our model are irreversible, and mean translocation rate J increases with both external force F and concentrations 
[EF-G], [EF-T], the mean dwell time T of ribosome in one translocation cycle decreases with both F and [EF-G], 
[EF-T], see Fig. H 



IV. DISCUSSION 

In this study, a model for ribosome translocation along messenger RNA is presented. In which, ribosome is assumed 
to translocate by two processes, active process and passive process. In each of the two processes, two sub-processes 
are included, translocation factor binding sub-process in which all the needed translocation factor are attached, and 
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mechanochemical sub-process in which one forward translocation step is made. The main difference between active 

process and passive process is that the transition rates of the mec/ianoc/iemicaZ sub-processes are different, which is due 
to the difference of free energy barrier. By this model, we found that, the mean translocation rate of ribosome increases 
with both the concentrations of translocation factors and the external force. Under high force, ribosome translocates 
mainly through active process, and at saturating concentrations of translocation factors, the translocation is limited 
by mechanochemical sub-processes. With the increase of external force, the mean translocation rate of ribosome 
increases from its lower limit to its upper limit. However, both of these two limits increase with concentrations of 
translocation factors. 
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TABLE I: Values of model parameters used in Figs. I4I7I which are obtained by fitting formulation ^ of mean translocation 
rate of ribosome along mRNA to the experimental data obtained in |23l] (see Fig. |4|. See also Fig. [2ja) and Eqs.fl}, (|9}, (|lip 
for definitions of the model parameters. 



Parameter 


Value 


Parameter 


Value 


kii 


13.78 (s^^ -fiM-^) 


ki2 


106.85 (s-^ -^iM"^) 


k2 


0.44 (codons-s^^) 


hi 


0.26 (codons-s"^) 




8.66 X 10~^ (s^^) 


kc 


41.28 (s^^) 


Fo 


0.66 (pN) 
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FIG. 1: Schematic depiction of the general model to describe helicase translocation along nucleic acid (a) and its corresponding 
two-line Markov chain (b). There are two possible processes to complete one step, active process (top, 1 — > 2 — >■ 1) and passive 
process (down, 1' — >■ 3 — >■ l'), with different translocation factor binding rates and mechanochemical transition rates ki, k2 and 
fei,fe3, respectively. These two processes can switch to each other with external force and state dependent rates, ko, kc, k'^, k'^. 
The needed translocation factors are attached during sub-process 1 — >■ 2 or l' — >• 3, and one forward step is completed during 
sub-process 2 — )• 1 or 3 — > l'. 




FIG. 2; Simplifications of the model depicted in Fig. [T] (a) The switch rates between active process and passive process are 
assumed to be state independent, k'^ = ko and k'^. — kc, and the rates of translocation factor binding sub-process in both active 
process and passive process are also assumed to be the same, k'l — ki. (b) states 1 and l' are combined to one state, and 
k'l — ki is assumed. Since the parameter numbers of the two simplified models are the same, but more information can be 
obtained from the one depicted in (a), in this study only the simplified model (a) is analyzed. 




FIG. 3: (a) The extension-force relation of DNA obtained by Tinoco and Bustamante in [43. (b) The free energy-force relation 
obtained by AGf = Jf^ z{F)dF (thick dotted line) with 2(F) the total extension under force F as plotted in (a). For F < 60, 
AGf can be well approximated by 16-F — 8 (thin solid line). 



11 



JO 

o 

T3 
O 
O 



B 
■5 
i_ 

c 
o 

'■I— • 
CC 

o 
o 

03 
CO 



c 

CO 

CD 




0.2 0.4 0.6 0.8 1 
[EF-G] (^M) 




I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I 




(b) 



■ ■■■■■■■■■■■■■■■ 



0.2 0.4 0.6 0.8 1 

[EF-G] (^M) 



i 



T 



...i 

(d) 



0.2 0.4 0.6 O.i 

[EF-T] (^M) 



FIG. 4; Theoretical results of ribosome translocation along hpValccso rnRNA [see Eqs. (O © ^-nd pip] with parameter 
values listed in Tab. [D The data points are experimentally obtained by Qu, et al 23]. In (b) (c) and (d), the solid curves are 
mean translocation rate J under high external force {F — 20 pN is used in the calculations), and dashed curves are values for 
low external force {F — 2 pN is used in the calculations). In (a) [EF-G] = [EF-T] = 1 fiM is used. 
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FIG. 5: Mean translocation rate J of ribosome along lipValccso bmRNA (see [23[ |) as functions of external force F, concen- 
trations [EF-G] and [EF-T]. Parameter values used in the calculations are as follows, (a) [EF-T] = 1 jiM, and reading from 
the bottom up [EF-G]=0.05, 0.1 0.2, 0.3, 0.5, and 1 fiM respectively, (b) [EF-G] = 1 ^M, and reading from the bottom up 
[EF-T]=0.01, 0.02, 0.03, 0.05, 0.1, and 0.6 fiM respectively, (c) [EF-T] = 1 fiM and (d) [EF-G] = 1 ^M. In (c) (d), reading from 
the bottom up, external force F =11, 13.5, 14, 14.5, 15, and 16 pN respectively. Other parameter values are listed in Tab. |T] 
The calculations indicate that the mean translocation rate J increases with both external force F and concentrations [EF-G], 
[EF-T]. 
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FIG. 6; Probability p\ + pi of ribosome at state 1 or l' [see Fig. Ofa)] as function of external force F (a) and concentration 
[EF-G] (c). (b) (d) are for probability pi +P2 of ribosome at state 2 or 3. The parameter values used in the calculations are as 
follows, (a) [EF-T] = 1 ^M, and reading from top downwards [EF-G]=0.01, 0.02 0.05, 0.1, 0.2 and 1 pM respectively. (b)[EF- 
G] = [EF-T] = 1 ^M, the curves are independent of [EF-G] and [EF-T]. (c) (d) [EF-T] = 1 fiM, and reading from the bottom up 
F =11, 13.5, 14, 14.5, 15, 16 pN respectively. For other parameter values, see Tab. ID 
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FIG. 7: Mean dwell time of ribosome in one translocation cycle as the change of external force [see (a) and (b)], concentrations 
of EF-G»GTP [see (c)] and EF-Tu»GTP»aa-tRNA [see (d)]. In addition to the parameter values listed in Tab. [H others are 
as follows, (a) [EF-T] = 1 fiM, and reading from top downwards [EF-G]=0.05, 0.1, 0.2, 0.3, 0.5 and 1 fiM respectively, (b) 
[EF-G] = 1 /^M, and reading from top downwards [EF-T]=0.01, 0.02, 0.03, 0.05, 0.1, and 0.6 fiM respectively, (c) [EF-T] = 1 fiM, 
and (d) [EF-D] = 1 /iM. In both (c) and (d), reading from top downwards F —11, 13.5, 14, 14.5, 15, 16 pN respectively. 



